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Abstract 



We discuss the non-linear growth of the excess kurtosis parameter of the smoothed density 
fluctuation field 5, S4 = [{S^) — 3(5^)^]/(5^)^ in an Einstein-de Sitter universe. We assume 
Gaussian primordial density fluctuations with scale-free power spectrum P{k) oc /c" and analyze 
the dependence of on primordial spectral index n, after smoothing with a Gaussian filter. 
As already known for the skewness ratio ^'3, the kurtosis parameter is a decreasing function of 
n, both in exact perturbative theory and in the Zel'dovich approximation. The parameter 5*4 
provides a powerful statistics to test different cosmological scenarios. 

Subject headings: Galaxies: clustering - large-scale structure of the Universe 

1 Introduction 

The study of the statistical distribution of the matter in the universe may be a way to address 
fundamental issues such as the origin and the formation of structures on large scales. 

The simplest and most usually accepted hypothesis, supported by the inflationary model, 
is that the very early distribution is Gaussian. In such a case, the connected A^-point corre- 
lations (and the A^-th order connected moments) with N > 2 are zero. However, even if the 
primordial fluctuations d are Gaussian, the non-linear time evolution will ensure that the mass 
density fluctuations become highly non-Gaussian (Peebles 1980; Fry 1984). It is important 
to understand the nature of the higher moments of the mass density induced by gravity in 
order to distinguish their effects from those of possible primordial non-Gaussian fluctuations: 
late-time phase transitions, cosmic string models and global textures are indeed models whose 
statistics may not be described by a Gaussian distribution (see e.g. Vilenkin 1985; Turok 1989; 
Scherrer & Bertschinger 1991). Moreover, variations of the inflationary model which lead to 
non-Gaussian primordial fluctuations have been recently discussed (see e.g. Salopek 1992). 

A powerful method to distinguish if the non-Gaussian nature of the matter distribution is 
intrinsic, is to analyze the growth of higher moments, like the skewness or the kurtosis (third 
and fourth connected moments, respectively) of the density fluctuation field. Peebles (1980) 
first showed that gravity induces, for an unsmoothed initial Gaussian density field, a skewness 
ratio Ss = {5^) / — 34/7, for any primordial power spectrum. Juszkiewicz, Bouchet, & 
Colombi (1993) find that the filtering operation actually introduces a dependence of ^'3 on the 
primordial spectral index n; in particular for scale-free spectra, 5*3 = 34/7— (n+3) for a top-hat 
window function, and a decreasing trend with n is recovered also for a Gaussian filter. Coles et 
al. (1993) investigate, using A^-body simulations, the growth of the skewness ratio to test the 
hypothesis of Gaussian primordial density fiuctuations against possible alternatives. Analytical 
expressions for 6*3 arising from non-linear evolution in perturbation theory have been worked 
out for arbitrary non-Gaussian models by Pry & Scherrer (1993) and Catelan & Moscardini 
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(1993). Finally, higher order moments are known to depend extremely weakly on the density 
parameter Q (see Martel Sz Preudling 1991; Bouchet et al. 1992; Bernardeu 1992). 

Here we analyse the dependence of the induced-by-gravity excess kurtosis of the density 
field, smoothed with a Gaussian filter, namely the parameter 5*4 = [{S^) — 3(5^)^]/ (5^)^, on an 
initial (scale-free) power spectrum P{k). To do this, we take advantage of the exact perturbative 
technique (Fry 1984; Goroff ct al. 1986) and the Zel'dovich approximation (sec Grinstein & Wise 
1987). The kurtosis describes features such as sharpness or strctchincss of the mass distribution 
and the extent of its rare-event tail. Moreover, it is possibly related to the initial sign of the 
skewness - as predicted in some non-Gaussian models (see e.g. Luo & Schramm 1993) - which 
is important for the final galaxy clustering pattern (Moscardini et al. 1991; Messina et al. 1992; 
Weinberg & Cole 1992). 

The layout of this Letter is as follows: in Section 2 we review the exact perturbative theory 
and the Zel'dovich approximation; in Section 3 we discuss the induced-by-gravity kurtosis 
parameter 6*4 of an initial Gaussian density field and its dependence on the primordial spectral 
index; we state the main conclusions in Section 4. 

2 Non— Linear Time Evolution 

We assume that present-day structures in the universe formed by gravitational instability from 
Gaussian fluctuations (5 in a pressureless fluid with matter density p — p{,[l + 5], where is the 
background mean density. The density fluctuation fleld 5 may be written as a Fourier integral, 
5(x, i) = (27r)~^/ dk5(k, t) e*'^'^ , where x and k are the comoving Eulerian coordinate and 

wavevector respectively, t is the cosmic time. The power spectrum P{k), defined as the 2-point 
correlation function in Fourier space (at a given time), (5(ki) 5(k2)) = {211)^ 6 0(^-1 + ^2) -P(^i) 5 
fully determines the statistics of the primordial Gaussian density field, whose variance is given 
by the expression (x^ = (1 /2Tr^) dk k^ P{k). 

To bridge the gap between observational data and theory, it is necessary to filter the fluc- 
tuation fleld 5 by means of a window function Wr , (5k(x, t) = J dy 6{y, t) H^Rdx — y |) . In the 
following we will adopt a Gaussian window function. The mass variance on scale R, a% , is 
related to the primordial spectrum by crj = (l/27r^) dk P{k) [WR{k)]^ , where WR{k) is 
the Fourier transform of Wr{x) . 

2.1 Equations of Motion: Perturbative Theory 

The time evolution equations for the matter density fluctuation 5(x, t) and the peculiar velocity 
fleld v(x, i) are the Poisson equation, the Euler equation and the continuity equation, i.e. 



(1) 



5oV + - (v V) V + - V = — V$ , 



a a a 



(2) 
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+ - V • (1 + (5) V . (3) 

a 

Here do = d/dt and spatial derivatives are with respect to x. We analyze these equations 
assuming an Einstein-de Sitter universe with vanishing cosmological constant. In such a model, 
the scale factor a is proportional to t during the matter dominated epoch, and the adiabatic 
expansion implies that Gtt Gpb — 1 . The quantity $ is the Newtonian gravitational potential. 
Combining the divergence of the Euler equation with the continuity equation, a second order 
differential equation for the density contrast 5 may be introduced (Peebles 1980) 



dtS + 2-doS-A7rGpbS^A7rGpbS' + — djda^ + dada(l + 5)v'^v'' . (4) 

The first term in the r.h.s. of Eq.(4) corresponds to the spherical collapse of an isolated proto- 
object, while the "geometrical" term in square brackets describes tidal and shear effects. The 
linear approximation, adequate when 5^ <^ 1, corresponds to dropping the r.h.s. of Eq.(4). The 
first order solution has the well-known self-similar form (considering only the growing mode) 

<5W(x,i) = D(i)5i(x), (5) 

where D{t) oc a{t) is the time growth factor of the mass fiuctuations. Higher order approxima- 
tions of the solution of Eq. (4) may be recovered if one expands the mass density fiuctuation field 
5(x, t) about the background solution S — (and v = 0), namely S — J2n ^i^^ — ^i^i): 
then solving the differential equation for any 5^"^ (Peebles 1980; Pry 1984). The perturbative 
expansion for 5 is (e.g. Goroff et al. 1986) 

5(x,t) = f;[D(t)r5„(x). (6) 

n=l 

The first term of the expansion corresponds to the linear approximation. Second and third 
order solutions have been discussed by Peebles (1980) and Fry (1984) respectively. We see that 
the scale factor D{t) acts as a couphng constant in this perturbative approach, since 

Here we review the exact perturbative technique to solve approximately the equations of 
motion of a pressureless gravitational fiuid up to third order in the density fiuctuation field. In 
particular, we will use the third order solution to compute the fourth order moment (namely 
the kurtosis parameter 5*4). We use exphcitly the same notation of Pry (1984). 

(i) Second Order Density Solution 

The second order contribution 5^"^^ is related to the second order gravitational potential ^^"^^ by 
the Poisson equation V^^'-^-* = AnGa'^phd^'^^ ; it is solution of the differential equation 




5(1) 2 + 
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+ 



(7) 



where A*^^) is the rescaled hnear gravitational potential defined by A*^^) = ^'^^'^ / AnGphO^ ■, for 
which V^A^-*^) = 5^^\ Each side of Eq.(7) is homogeneous in powers of t. For an Einstein-de 
Sitter universe, noting that D"^ oc 5^^^ = D'^{t) 52(x), the solution of Eq.(7) is given by (Peebles 
1980, §18) 

5^'^ = \ + 9. A« + ^ 9«9^A« 9.9, A(^) . (8) 

Note that (5(2)) = o, i.e. mass is conserved to second order. We see from Eq.(8) that the 
behaviour of 5 at second order is non-local: the mass fiuctuation at the position x depends on 
initial perturbations at other positions via A^^^ . Physically this means that, unlike the linear 
local case, when density fluctuations grow in amplitude their spatial dependence, in comoving 
coordinates, changes. Thus, the gravitational field changes direction and particles are not 
accelerated in a fixed direction, as it occurs in linear regime. The last term in right-hand side 
of Eq.(8) corresponds to the velocity shear contribution. 

We can obtain the Fourier transform ^^^-^ directly from the differential equation (8). Explic- 
itly (time dependence is understood). 



^(^)(k) =/7^ J(^)(k',k-k') 5«(k')5«(k-k') , 



(27r) 



(9) 



where we have defined the kernel 



kk^^ 

kk' 



(10) 



(a) Third Order Density Solution 
The third order approximation 5^^^ is solution of the differential equation 

+ {dj^'^ 9,$(2) + dj^'^ + a-' d^d^ (2 '^v^^'^^ + 5^'^ v^'^'^ v^^^^) . (11) 

In this case, working directly in Fourier space is much simpler. Since 5'^'^^ oc oc t^, we have 
^2^(3) ^ 2da"^9o5(^) - A-nGpi^l^^^ = Ad^^^t^ ; by using the second order results it is not difficult 
to show that (Fry 1984) 



5(^)(k) 



rfki d'k2 d\ii 



(27r)6 

where the third order kernel is 

J(3)(ki,k2,k3) = j(2Hk2,k3) 



(5^(ki + k2 + k3-k) J(3)(ki,k2,k3) 5«(ki)^W(k2)^('Hk3) , (12) 



11 ki • (k2 + k3) 4 kJCi k • (k2 + k3) 
3^3 (k2 + k3)2 ^9 (k2 + k3)2 



+ 
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2 k • ki k ■ (k2 + ka) (k2 + kg) • kg 1 k • ka k ■ kg 

9 (k2 + k3)2 ki ^ 9 ' ^ 

It is clear from Eq.(ll) that, in order to derive the solution 6'^^\ we need to know explicitly the 
second order velocity v^^^. This is given by 



(2) D 



5(1) VA« -2VA(2)] , (14) 



where A^^^ = ^^'^yAnGpb is the rescaled second order gravitational potential. [In the solution 
(14) we neglect an additive homogeneous term whose divergence is zero.] In an Einstein-de 
Sitter universe v^^^ ~ t, slower than 5^^^ ~ t^/^. We stress the fact that v^^^ is not parallel to 
the second order acceleration [oc — VA*^^)]: this is a consequence of non-locahty. Finally, we 
note that v^^^ is known only once 5^"^^ is known. 

(in) General Solution 

As Goroff ct al. (1986) have shown, the (Fourier transformed) n-th order mass fluctuation term 
may be represented in integral form as 

ny (^^i(k,) I [(27r)%(5:k,-k)] jW(k„...,k„). (i5) 

The presence of the Dirac delta function comes from momentum conservation in Fourier space. 
The kernels J^"^ are symmetric homogeneous (with degree 0) functions of the wavevectors 
ki, . . . , kjj , describing the effects of non-linear collapse (tidal and shear forces). In general J*^"^ 
arc very complicated for n > 3. [A discussion of the properties of the kernels J*-"-* is given in 
Wise (1988). Explicit recursion relations with their Feynman diagrammatic representation are 
given by Goroff et al. (1986) and Wise (1988).] 



2.2 Zel'dovich Approximation 

In the Zel'dovich approximation (Zel'dovich 1970) the motion of particles from the initial co- 
moving (Lagrangian) positions q is approximated by straight paths. The Eulerian position at 
time t is then given by the uniformly accelerated motion 

x[q,i]=q + L>(i)S(q), (16) 

where D{t) is the growth factor of linear density perturbations and S(q) is the displacement 
vector related to the primordial velocity field. The Zel'dovich approximation provides an exact 
solution of the equations of motion for one-dimensional perturbations, and reduces to the 
linear approximation when 5 and v are small. In general, the Zel'dovich approximation is 
not an exact solution of the equations of motion, in that in a finite time particles converge 
into singular regions of infinite density (caustics), and the map in Eq.(16) becomes multi- 
valued. The treatment of these formal singularities is the main problem to solve during the 
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highly non-hnear stage of structure formation. Smoothing on a suitable scale R partially 
solves this problem. Grinstein & Wise (1987) give an Eulerian representation of the Zel'dovich 
approximation by a diagrammatic pcrturbative approach similar to that of the previous section. 
They showed that the n-th order pcrturbative corrections 5„(x), when the density fluctuation 
fleld 6 is evolved according to the Zel'dovich approximation, are such that 

°° ^8 (9 r 1 

^(x, ^) = E ^ mr EgT-'-'gr- [S^. ■ ■ ■ Sh.] ■ (17) 

n=l "'■ [hn]=l 

Here J2[hn] = ' ' '^hn- Note that the first term recovers the linear approximation, in that 
S = vi, where Si{x.) = — V-Vi . The pcrturbative expansion for S of Eq.(17) simply corresponds 
to calculate different symmetric kernels J^a- These can be written in the following compact 
form 

jak„...,k„) = in^, (18) 
^- h=i % 

where k = Y.h=i ■ Note that, unlike the pcrturbative case, the kernels J^^ are manifestly 
symmetric by construction. 



3 Kurtosis of the Density Field 

In this section, we compute the induced-by-gravity kurtosis parameter S4 of an initial Gaussian 
density field in a fiat universe. The lowest order non-zero contribution to 5*4 is 

{S^'^YS, = 6 5(^)2) +4(5«35(^)). (19) 

From Eqs.(lO) and (13), after tedious but straightforward algebra, one finally gets the integral 
expression of the kurtosis ratio. 



24 rdkidkidks 
~^ J (2^09 



Pih)P{h) [P(A;3) J(=^)(ki,k2,k3) 



+ 



+2 j(^H-k2, k2 + kg) J(^)(ki, k2 + kg) P(|k2 + kgl)] . (20) 

The expression in the Zel'dovich approximation is obtained by substituting the kernels J^'^^ 
with the corresponding J^]^ in the previous relation. The angular integrations in Eq.(20) may 
be analytically performed. One obtains S4 — 60,712/1,323 ^ 45.9 in the pcrturbative case 
(Pry 1984; Bernardeu 1992) and ,^4 = 88/3 ^ 29.3 in the Zel'dovich approximation. These 
unsmoothed-case results are independent of the primordial spectral index, as is the skewness 
ratio 5*3. 

The kurtosis of the smoothed density field Sr in the exact pcrturbative case is 
S^iR) = ^ / '^^' f'^J^' Wnih) Wnih) Wnih) WnHk, + k2 + kgl) x 
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P{h)P(k2) [Piks) J^^\k,MM) + 2j(2)(-k2,k2 + k3) j(2)(ki,k2+k3)P(|k2 + k3|)] ; 

(21) 

again the substitution J*^") J^*^ leads to the corresponding expression in the Zcl'dovich 
approximation. Note that P{k) completely describes the process of growth of mass fluctuations 
from Gaussian initial perturbations. 

4 Discussion and Conclusions 

We calculate the previous integrals, by an Adaptive Multidimensional Monte Carlo Integration 
subroutine, for scale-free power spectra P{k) oc A;" with n in the range — 3 < n < 1. Due to the 
assumed primordial scale-invariance, 5*4 only depends on the primordial spectral index n, and 
not on the scale R. In Figure 1, we plot the kurtosis ratio S4 of the Gaussian-smoothed density 
field versus the spectral index n, for both the perturbative and Zel'dovich approximations. 

It clearly appears that S4 strongly depends on the primordial spectral index n. In particular 
the kurtosis parameter is a decreasing function of n . This is also confirmed by Bernardeau 
(1993), who, applying the exact perturbative technique, finds a similar trend in the case of 
top-hat filtering. Anticorrelation with the amount of small-scale power has also been found 
for the skewness parameter S3 (Fry 1984; Juszkiewicz, Bouchet, & Colombi 1993), and it seems 
therefore a general property of higher order moments of the smoothed matter distribution: 
larger values of n correspond to higher power on small scales, where the filtering operation 
acts. 

Furthermore, the Zel'dovich approximation underestimates the induced-by-gravity S4 w.r.t. 
the rigorous perturbative one, as already noted by Grinstein & Wise 1987 (although in the 
framework of the "standard" cold dark matter model). This is hardly surprising, since the 
Zel'dovich approximation fails to fully describe the gravitational effects causing particle tra- 
jectories to depart from their original directions. It results that the Zcl'dovich approximation 
makes higher orders in the perturbative series for 6 smaller on large scales than they actually 
are (see e.g. Wise 1988). 

We stress that, due to the importance of the smoothing procedure, but also to its arbitrari- 
ness, one must be cautious about making quantitative comparisons between our results and 
both observational and A'^-body data. A Gaussian filter is used for instance by Saunders et 
al. (1991), who however estimated only the skewness in the QDOT-IRAS catalog. The only 
observational estimates up-to-now of the kurtosis of the galaxy distribution are obtained by 
counts-in-cells, corresponding to a top-hat filtering (Gaztanaga 1992; Bouchet et al. 1993; Fry 
& Gaztanaga 1993); a similar method is applied by Lahav et al. (1993) and Lucchin et al. (1993) 
in A^-body simulations with scale-free power spectra. Our results can be used to constrain both 
the probability distribution and the power spectrum of primordial density fluctuations and thus 
help to select the most reliable mechanism for the generation of large-scale structures. 
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Figure caption 

Figure 1. The kurtosis ratio 6*4 of tlie density field for power-law spectra P{k) oc /c" and 
Gaussian filter versus the primordial spectral index n, for both the perturbative (squares and 
solid line) and Zel'dovich approximations (triangles and dotted fine). Note that the values at 
n — —3 correspond to the unsmoothed cases. Error bars refer to the associated uncertainty 
estimate from the Monte Carlo Integration. 
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